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Abstract. Computing modular coincidences can show whether a given substitution system, 
which is supported on a point lattice in R d , consists of model sets or not. We prove the 
computatibility of this problem and determine an upper bound for the number of iterations 
needed. The main tool is a simple algorithm for computing modular coincidences, which 
is essentially a generalization of Dekking coincidence to more than one dimension, and the 
proof of equivalence of this generalized Dekking coincidence and modular coincidence. As a 
consequence, we also obtain some conditions for the existence of modular coincidences. In a 
separate section, and throughout the article, a number of examples are given. 



1. Introduction 

Substitution tilings, or substitution point sets, are examples of objects which show high 
order at both short and long range, even if they are not periodic. This manifests itself in 
the diffraction spectrum of such a tiling. Like crystals or periodic tilings, some nonperiodic 
substitution tilings show a diffraction pattern consisting of bright spots — the "Bragg peaks" 

— only. Such tilings are called pure point diffractive. In general, it is hard to decide whether 
a given substitution tiling is pure point diffractive. For a special class of tilings the question 
can be answered for any given example in principle: If the set of difference vectors between 
tiles spans a discrete point lattice, then there is a condition which is equivalent to pure 
point diffractivity ( |LMS2j . see also Theorem Yd. 31 below). This class covers many well-known 
examples, such as the chair tiling and the sphinx tiling (being pure point, see, e.g., |Soj 

BMSj, [E], |LM| ). or the Thue-Morse sequence and the table tiling (being not pure point, 
see, e.g., |Ma| . |K] . |R]). In fact, the condition of |LMS2j may be difficult to check in some 
cases. Here, we devise a simple algorithm and show the equivalence between the output of 
this algorithm and the condition mentioned. As a consequence, we obtain results about the 
computability of the problem and some conditions for the existence (resp. nonexistence) of 
modular coincidences. 

The paper is organized as follows: Section [2 contains notations and basic terms used in this 
text, the definition of lattice substitution systems and related facts and definitions. Section 
13 is dedicated to the main result of |LMS2j ( Theorem 13. 3j) . together with the definitions of 
model sets and modular coincidence. Since we used the notation of LMS2j wherever it is 
possible, readers who are familiar with |LMS2j may skip Sections El and 01 In Section |IJ we 
give the construction of the coincidence graph, which is a generalization of Dekking coinci- 
dence [D], and establish the equivalence between this coincidence and modular coincidence 
(Theorem I4.4|) . From this result, we obtain results about the computability of the problem 
(Theorems 14.51 f4. 12j) . as well as some conditions (Corollaries 14. f>| I4.7| I4.9() for the existence 
(resp. nonexistence) of modular coincidences. Section is a brief remark about the best 
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upper bound for the costs of computing modular coincidences. Section is a collection of 
examples which illustrate terms and statements of preceding sections and show the usefulness 
of our results. Finally, three appendices show the connections between this text and Dekking 
coincidence, the overlap algorithm in |Soj and automatic sequences, respectively. 



2. Basic Definitions, Substitution Systems 

Let us fix some notation. By N we denote the set of all positive integers. The sum of two sets is 
to be understood point- wise, i.e., M + N = {n + m\m € M, n € N}. The closure of a set M is 
denoted by cl(M), the interior of M by int(M). The cardinality of A is denoted by j^A (where 
#(0) = 0). By M d we denote the closed unit ball {x | < 1} in M. d . A tile is a nonempty 
compact set T C M. d with the property that cl(int(T)) = T. A tiling of 13L d is a collection of 
tiles {T n } ngN which covers M. d and contains no overlapping tiles, i.e., int(T fc ) n int(T n ) = 
for k 7^ n. A tiling T is called nonperiodic, if the only solution of T + x = T is x = 0. If 
the tiles of a tiling belong to finitely many translation classes [TJ, [T 2 ], . . . , [T m ], then the 
class representatives T X ,T 2 , . ■ ■ ,T m are called prototiles. If necessary, one may consider two 
congruent prototiles as different by labelling them with their type. Given a set of prototiles 
Tj, one can describe tilings — or parts of tilings — by a set 

(1) {T 1 + t, « , T x + 4 1} , . . . , T 2 + vf ] , T 2 + 4 2) , . . . , T m + , . . . } € M d ) , 

or shortly by {T x + T 2 + . . . , T m + V m } — where V i contains all vfi of — or even 
shorter by a position vector V = (V 1 , V 2 , ■ ■ ■ , V m ) T . 

Nonperiodic tilings of long range order can easily be generated by a substitution. The idea 
is, basically, to give a set of prototiles and a rule how to enlarge them, and then to dissect 
them into tiles congruent to the original prototiles (cf. Example 16. 7|) . Formally, this can be 
done by a matrix function system, resp. by a substitution system. 

Definition 2.1. A matrix function system (MFS) is an m x m-matrix $ = ($y)i<*J<m; 
where each is a finite (possibly empty) set of mappings x t— > + o ijfc; where Q is an 
expansive linear map, 1 < i,j < m, 1 < k < a ijfc € 

j4n MFS is called primitive if its substitution matrix = (#(^f?))i<i j<m (a^ x wi- 
matrix with non-negative integer entries) is primitive, i.e., S k has only strictly positive entries 
for some k E N. 



Note that we can define Q k (for k £ N) inductively by the composition of maps, i.e., 
= ^^ = (^° $ )i<i,i<m = (Uti^k ° (*)«)i< i , i <m where 



if ($ fc ) i£ = or $ ej = 0, 

{5 / I 9 e ($ fc )ii, / £ *,.■} otherwise. 



Clearly, we have S($ k ) < «S'( < I ) ) fc . Note also that the linear part of the affine mappings in <£ fe 
is given by x 1— > Now we can apply the MFS $ to a position vector V. I.e., if V = 

(Vi,..., V m ) is a position vector as above, then ®V = (W\, ... , W m ), Wi = \JjLi = 

Definition 2.2. A substitution system is a pair (V, $), where V = (V 1 ,V 2 , . . . ,V m ) T is a 
fixed point of the matrix function system $ (i.e., = V), such that all V { are mutually 
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disjoint and the unions in 



m m 



Vi=\j w = u u m 



are pairwise disjoint for all 1 < i < m. 

If supp(V) = IJ™ x Vi is a lattice, then (V, <E>) is called lattice substitution system (LSS) with 
m components. 

Now the relation between substitution systems and substitution tilings becomes clear: A set 
of tiles is represented by a vector V whose entries Vi are point sets (where each point encodes 
the position of one tile). The substitution applied to this set of tiles yields a larger set of 
tiles, which is represented by the vector <&V. 

By the action of a tile substitution, each tile T is enlarged by Q, then dissected into some 
tiles of type [T x ], . . . , [T m ]. If T € [TJ, then QT is dissected into #( < I ) li ) tiles congruent to 
T 1; #(&2i) thes congruent to T 2 and so on. The translational part of the affine maps in & 
determines the exact location of these tiles inside the enlarged tile. In this way, (T, <I>) can be 
considered a substitution system. Wherever we will emphasize the fact that (T, $) describes 
a tiling, we will call T a substitution tiling and $ a tile substitution. For examples of the 
relation between tile substitutions and substitutions on point sets, we refer to Example 16.71 
The next lemma reads essentially: The translate of an LSS (V, <£) by a t £ supp( V) is again 
an LSS. 

Lemma 2.3. Let (V, <£) be an LSS, t G supp(V), and let (V - t,& + Qt - t) be defined by 
V-t = (Vt-t, . . . ,V m -t) and^+Qt-t = ($ ii +Qt-%, ^ij+Qt-t = {Lp+Qt-t\ip € 
Then (V — t, <& + Qt — t) is also an LSS. 

Proof. A simple calculation (cf. Def. E2J) yields: ($ + Qt - t)(V - t) = V - t if and only if 
&(V) = V. Since t £ supp(V), it follows that supp(V" — t) is a lattice, thus all conditions for 
lattice substitution systems are satisfied. □ 

If (V,&) is an LSS, then the set V in Definition 12.21 consists of point sets Vi. Rather than 
considering V as a subset of M. dm , we may think of V as a point set in M^, where each point 
has got one color (out of m colors altogether). Such a single point of V with color j, located 
at x £ M. d , is of the form 



which is a somewhat inconvenient notation. Therefore, we will avoid to consider such single 
points as much as possible. If such a single colored point located at x of type j is considered, 
we use the notation 'V n {x}, x is of type j' (or of color j). 

In this paper we consider only primitive LSS (V,<&). This implies that each V i is a Delone 
set, and we call V a multi- component Delone set (see |LMS1| Section 2], where the name 
Delone multiset is used; however, the name "multiset" in this context is misleading). 

Definition 2.4. A set V C M. d is a Delone set (or Delaunay set) if it satisfies the following 
two conditions: 



{(0,...,0, {x} ,0,...,0)} 



jth entry 
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(1) V is uniformly discrete: There is a value r such that each open ball of radius r contains 
at most one point ofV. 

(2) V is relatively dense: There is a value R such that each closed ball of radius R contains 
at least one point ofV. 

Note that the maximal such radius r is the packing radius for V, while the minimal value of 
R is the covering radius of V. 

We should mention that every primitive substitution tiling gives rise to a primitive substitu- 
tion system in a quite natural way, see Q and comments there. In the other direction, this 
is less obvious. In the case of LSS, one simple way is the following: Assign to each point x 
of type i in a multi-component Delone set a tile x + F of type i, where F is a fundamental 
domain of the lattice. Then it is trivial to read off the substitution for the tiling from the 
MFS for the multi-component Delone set. Note that we do not require the tile-substitution 
to be self-similar in the sense of |Soj . For details about the subtleties of the relation between 
point-substitutions and tile-substitutions, see |LMS2| Section 3.2]. The property which de- 
scribes this relation is given in the following definition. It works both for tilings and for 
multi-component Delone sets. In order to express it, we denote by C r (x, A) the finite portion 
of the structure A in a ball of radius r around x, the r -cluster at x. If A is a tiling, then 
C r (x,A) := {T G A\Tn(rM d + x) + 0}. If A = (At, . . . , A m ) is a multi-component Delone 
set, then C r (x,A) := (A 1 n (rM d + x), ...,A m n (rM d + x)). 

Definition 2.5. Let A,B be two multi-component Delone sets, or two tilings, or one multi- 
component Delone set and one tiling. B is locally derivable (with radius r) from A, if, for 
all x, y in M. d , one has: 

C r (x,A) = C r (y,A) + (x-y) =^ C (x, B) = C (y, B) + (x - y). 

If B is locally derivable from A and A is locally derivable from B, then A and B are mutually 
locally derivable (MLD), see (HSU- 

For examples of a tiling and its corresponding substitution system being MLD see Example 
16.71 Another property which is frequently used is the following. 

Definition 2.6. A tiling A, or a multi-component Delone multiset A, is called repetitive, 
if for every r > there is R > such that each R-cluster in A contains a translate of each 
r -cluster in A. 

3. Model Sets 

The objects we consider in this paper are primitive lattice substitution systems (V,&) and 
tilings which are MLD to a primitive LSS. To exclude certain pathological cases, we assume 
further that all r-clusters in V are contained in some & k (v), where k € N and v is a single 
point of some type i. In particular, this ensures that V is repetitive. Moreover, the dynamical 
system associated to V is always minimal, see the remark at the end of this section or |Soj 
for details. 

A central question is whether such a substitution system consists of model sets or not; and 
whether the diffraction spectrum of such a substitution pattern is pure point or not. If V 
is fully periodic (i.e., there are d linearly independent translation vectors such that V is 
invariant under these translations), this is trivial. The question becomes interesting if V is 
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not fully periodic. A theorem by Schlottmann |Scl Theorem 4.5] states that all regular 
model sets are pure point diffractive. We first recall what a model set is, see Mo_. 

Definition 3.1. A cut and project scheme consists of a collection of spaces and mappings as 
follows: 

R d ^L. R d x H ^ H 
U U U dense 

L i— i L — ► L* 

where H is a locally compact Abelian group, L cM. d xH is a lattice, i.e., a discrete subgroup for 
which the quotient group (M. d x H) /L is compact, and tt, 7r int denote the canonical projections, 
where ir int (L) = L* is dense in H and the restriction ir\^ is injective (and therefore bijective 

as map from L to L = vr(L)J. We call R d the direct space (or physical space j and H the 
internal space. 

A model set (or cut and project set) in M. d is a subset ofM. d which (up to translation) is of 
the form A(Q) = {tt(x) \ x 6 L, 7r int (x) £ Cl} for some cut and project scheme, where O C H 
has non-empty interior and compact closure. We call the window of the model set. 
A model set A(Q) is regular if the boundary dQ = cl(f2) \ int(fi) is of Haar measure 0. 

Note that one can define the so-called star map x* of x € L by x* = n iat (ir\~£ (x)), justifying 
the notation L* above. 

To an LSS, one can naturally associate a cut and project scheme (see Definition 13.21 below 
for the definition of V). The internal space is then given by H = (L/L') x li m< _ fc L'/Q k L', 
where lim^ L'/Q k L' denotes the Q-adic completion of V (a profinite group), see [LMS2 , 
Section 5.1], and Examples 16.11 and 16.21 Furthermore, the lattice L also arises naturally as 
diagonal embedding of L in M. d x H, i.e., L = {(t,t) \ t G L}. So only the existence and 
shape of a window is in question. For this, a central role is played by the notion of a modular 
coincidence. 

Definition 3.2 (cf. [LMl ILMS2j 1. Let (V, $) be an LSS and L = supp(F). Let V := 

Li+L 2 -\ h L m , where Li := (V{ - Vi)z- 

For a E L, let 

^ij[a] := {if G $^ [ ip{y) = a mod QL', where V j C y + L'} 
= { (p e<Z> ij \<p(V j )ca + QL r }. 

Furthermore, let 

*H:= [J %[«]• 

l<i ,j<m 

(V, $) admits a modular coincidence relative to QL', if $>[a] is contained entirely in one row 
of <3? for some a € L, i.e., if (a + QL') C for some 1 < i < m. 

For an illustration of this concept, see Example 16.41 with an explicit calculation of a modular 
coincidence in a particular case. (At this point, consider just the first part of the example, 
before the first occurrence of ^f). 

Remark: Since Q{L) C L, one also has Q(L') C V (see |LMS21 Section 5.1]). Furthermore, 
L/L' is finite, i.e., L' partitions L in finitely many cosets. By the definition of L' , each V- is 
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a subset of one such coset (it is possible that more than one Vj are subsets of a fixed coset). 
Clearly, we have V- C y + L' for y G V-. Since (p G $ i • is affine, we get ^piV-) C y?(y + L') = 
tp{y) + QL' . Therefore, $[o] partitions $ into congruence classes induced by the cosets L/QL' . 
Note that 3>[a] ^ for all a G L. Note also that even in the one-dimensional case, not all 
have to be equal, cf. Example 16.11 

The answer to the questions about an LSS being a model set and being pure point diffractive 
is now given in the following theorem. 

Theorem 3.3. |LMS21 Theorem 5.12] Let (V,$) (with V = (Vx,...,V m )) be a primitive 
LSS, L' = L\ H + L m , where Li := {Vi — Vi)%. Then the following claims are equivalent. 

(1) Each Vi is a regular model set. 

(2) There is a k G N such that a modular coincidence relative to Q k L' occurs in <3? fe . 

(3) V is pure point diffractive (meaning that each is pure point diffractive). □ 

Essentially, this theorem tells us that it is sufficient to compute powers of until there is a 
map ip which occurs in one row of <& k only. Alternatively, it follows from |LMS2j that it is 
enough to find one sublattice a + Q k L' which lies entirely in a single V^, to ensure that each Vj 
is a regular model set. At this point, it is difficult to estimate an upper bound for k, or even 
to decide whether the derived algorithms will terminate, assuming that V is not a model set. 
An answer is given by Theorem 14.51 below. 

Note that the above construction does not yield the "minimal" possible internal space, see 
Example 16.21 

Remark: To a multi-component Delone set V (or a tiling T), as obtained via a primitive 
substitution, one can associate a compact metric space Hy (X^-). Here, H v := cl({V — t \ 
t G i.e., it is the closure of the set of all translates of V in the space of all multi- 

component Delone sets, where the metric yields the following definition of closeness: Two 
multi-component Delone sets V 1 ,V 2 are close to each other if they coincide, after a small 
shift, inside a large ball around the origin. This compact space together with the translation 
action M. d yields a topological dynamical system (X^,IR rf ). If the eigenfunctions of this action 
span a dense subspace of L 2 (K V ), then we say that V has a pure point dynamical spectrum. 
In the cases considered here, V (resp. T) has pure point dynamical dynamical spectrum iff 
it is pure point diffractive. This equivalence is often used in the literature. In fact, this 
equivalence plays an essential role in the proof of Theorem 13.31 For details on dynamical 
systems and this last connection see |LMS1| ILMS2| IU| IBLj . 

4. Main Results 

In connection with modular coincidence, we define the following subsets of {1, . . . , m): 

* fc [o] := {i\(a + Q k L') n V t + 0} (a G L/Q k L') 
The following properties are obvious: 

(1) %[a] C {1, . . . ,m} and 1 < < m. 

(2) a + Q k L' C (J ^ 

ie* fc [a] 
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(3) (V, $) admits a modular coincidence relative to Q k L' iff there is an a G L/Q k L' such 
that_#(tf fc [a]) = 1. 

(4) For 2 G L/L', denote by z an element of the associated coset. Then ^q[z] is a partition 
of {1, . . . , m} as z runs through L/L' . 

(5) By Q fe+1 L' C Q k L' we have * fc [a] = U *fc+iH- 

2=a mod Q k V 

Already these properties can sometimes be used directly to decide on the (non-)existence of 
modular coincidences, see Example 16.31 

Note that the sets ^? k [a] can be calculated iteratively: 
^ k+1 [b] = {i\(b + Q k+1 L')nV^0} 



DefJO 



I 



6 + Q fc+1 L')n[U |J (piVj)] + 

3 <P£® tj [b] 



U U {(b + Q^L'Jn^)) +0 

3 f^ijlb] 



U U P [[a + Q k L'j n VA ± 0, where tp{a) = b mod Q k+1 L' 

3 <pe*iM 



= {i | there exist j, ip G &ij[b] and a G L s.t. (/5(a) = b mod Q k+1 L' and j G \& fc [a]}, 

where the fourth equation is due to the following: If (6 + Q k+1 L') n y(Vl) 7^ 0, then there 
exists an a such that 93(a) = 6 mod Q k+1 L'. 

For a fixed VE^a], we call each which is determined by the previous equation a child of 

Conversely, for a fixed we call all VPjJcz] which appear in the previous equation 

the parents of ^ k+1 [b}. 

Example 16.41 shall illustrate this concept. This example also inspires the following notation: 
We fix a system of representatives (which includes 0) in L of L/QL, denoted by a i . Then we 
can code every a G L by a finite sequence a = (a 1; a 2 , ■ ■ ■ , a n ) (the Q-adic expansion of a), 
where a G a 1 + QL, a — a 1 G a 2 + Q 2 L, a — a 1 — Qa 2 G 03 + Q 3 L, etc. 

Definition 4.1. ^4 primitive LSS (V, $) is called admissible, i/ a// mappings in /iave 
translational parts of the form . . .), where is an element of a fixed system of 

representatives for L/QL. 

Note that, by the definition of an LSS, all such translational parts appear exactly once in each 
column of & in the admissible case. The admissible primitive LSS are those to be considered 
in the sequel. Later we shall discuss nonadmissible substitution systems (Theorem I4.12[) . For 
an example of a nonadmissible primitive LSS, cf. Example 16.51 

In the following, we introduce a simple method for computing modular coincidences in the 
admissible case. 

Definition 4.2. The coincidence graph G(v,$) °f a primitive LSS (V, <5) is a directed 
(multi-) graph, where the vertices represent sets ^ k [a\. If^ k [a] and ^[6] are equal as sets, 
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then both are represented by the same vertex. 

We mark the [L : L'\ vertices ^ [a] as "base points" (e.g., these vertices are represented by a 
double circle instead of an ordinary circle). We draw an edge from v\ to V2 iff V2 is a child 
of v\ . Therefore, there are at most j det Q \ outgoing edges at each vertex. 

Of course, one can also label the edges by the element z € L' /QL' which induces the par- 
ent/child relation. This is useful if one wants to analyze the substitution further, e.g., de- 
termine the windows of a model set explicitly. In this case, multiple edges are possible, and 
there are exactly | detQ| outgoing edges at each vertex. For examples of coincidence graphs, 
see Examples 16.71 16.91 

The practical construction of a coincidence graph works as follows. Let us introduce the 
following notation: For an (admissible) LSS (V,<1>) and fixed j, denote by &(j) z the letter i 
for which Qx + z S (where z € L/QL'). By definition, this is (for every j and z) exactly 
one letter. Then we obtain the coincidence graph as follows: We start with the base points 
^oN ( a ^ L/L'). Suppose {s 1 , . . . ,s £ } = ^ [a] for some a, then there is a directed edge to 
the vertex representing {^(s^, . . . &(sf) z } for each z. From these new vertices, we proceed 
further (to the vertices { < I )2 (s 1 ) ;Z , . . . & 2 (sg) z }) until no further edges and vertices have to be 
added. This graph is then the coincidence graph. 

Lemma 4.3. Let (V,$) be an admissible primitive substitution system. If 1 5 k [a] = \f^[c], 
then ^/Ja] and \^[c] have the same set of children. 

Proof. Case L = L'\ Here a in ^ k [a] is actually a representative of L/Q k L, wherefore we can 
write (a ai , . . . , a ak ) for a. Similarly, we can write (a Cl , . . . , a c ) for c. An (admissible) map 
ip(x) = Qx + (a ijfc , 0, 0, ... ) yields (with these identifications) <p(a) = (a ijfc , a ai , . . . , a ak ) and 
93(c) = (cKyj., a Cl , . . . , di Ci ). Now, since ^ k [a] = ^[c], the same maps operate on a + Q k L 
as well as on c + Q l L, and we get ^ k+1 [t + Qa] = [t + Qc] (where t is an associated 
translational part). Furthermore, we observe in the admissible case that each child has exactly 
one parent ( 1 ^ k+1 [(a 1 , . . . , a k+1 )] has parent ^ k [(a 2 , • • • , Ofc+i)])- This proves the lemma for 
L = V. ' 

Case L ^ V: Since each set ^ k [a] determines the representative of L/L' uniquely, we can 
simply index the position by the Q-adic expansion relative to L'. This reduces this case to 
the previous one. □ 

It follows that we need to compute children only as long as we find sets that did not occur 
before. This justifies the representation of equal children ^? k [a] = ^[ft] by a single vertex in 
Def. 14.21 As a consequence, we obtain the following result. 

Theorem 4.4. Let (V, $) be an admissible primitive LSS. Then, (V,&) admits a modular 
coincidence iff at least one vertex of the coincidence graph G^y,^) represents a singleton set. 
The length of the shortest path from one of the vertices ^q[o\ to this singleton set is the 
minimal k such that a modular coincidence relative to Q k L' occurs in $ fc . □ 

The claim follows immediately from Lemma 14.31 and the construction of ^f. This result yields 
an upper bound for the number of iterations needed to compute modular coincidences. 
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Theorem 4.5. Let (V, <&) be an admissible primitive LSS with m > 2 components. If <& fc 
admits no modular coincidence for k < 2 m — m — 2, then <3? fc admits no modular coincidence 
for any k G N. 

Proof. As described above, it is possible to determine the sets ^ k [a] iteratively. The LSS 
( V, 3> fc ) admits a modular coincidence iff there is an a such that Vl/jja] is a singleton set. 

From Lemma l4.31 it follows that we need to compute children until no new sets appear. Since 
every child is a nonempty subset of {1, 2, . . . , m}, after at most 2 m — 2 iterations all possible 
sets have occurred. (We count the determination of ^f [-] as the zeroth iteration step. So 
corresponds to $ .) Moreover, after at most 2 m — m — 2, iterations either a set with 
just one element occurred, or, by Lemma 14.31 it will never occur. By Theorem 14.41 this is 
equivalent to the claim of the theorem. □ 

As a consequence of Theorem 14.41 we obtain some conditions for (V, <£) admitting a modular 
coincidence. By Theorem 13.31 this translates to necessary conditions for V consisting of 
model sets, resp. being pure point diffractive. 

Corollary 4.6. Let (V, <&) be an admissible primitive LSS. If there is an a such that ^ [a] 
is a singleton set, then V consists of model sets. 

Proof. This is an immediate consequence of Theorem 14.41 \P [a] is already a singleton set, 
thus, by the construction of we obtain Vi = a + V for some 1 < % < m. By the remark 
following Theorem 13.31 V consists of model sets. □ 

Corollary 4.7. Let be an admissible primitive LSS, such that no * s a singleton 

set. If V consists of model sets, then $ admits at least one pairwise coincidence in some 
VPq[ci], i.e., there are j ^ k such that j, k G VPqM f or some a and <&ij n <&ik 7^ for some i. 

Proof. Assume that $ admits no pairwise coincidence in any ^o[ a ]- Therefore, and because 
any set ^ [ a ] contains more than one element, any child ty^Qa + z] contains more than one 
element. Inductively, no set VPjJ-] is a singleton set. □ 

Surprisingly, the last corollary applies to a large number of substitutions found in the litera- 
ture, showing that they do not consist of model sets, since they admit no pairwise coincidence 
(cf. Examples lfi.3( lfT7| Ifi.Hjl . In |Fra2| Def. 1.2], the author defines bijective substitutions. 
This definition reads in our notation as follows. 

Definition 4.8. Let (V, $) be an admissible primitive LSS. If each ip occurring in <I> appears 
exactly once in each row, then (V, $) is called bijective substitution (system). 

The following statement is mentioned in |Fra21 Section 3.3], but without proof. Here we 
obtain the result quite easily. 

Corollary 4.9. Let (V, be an admissible primitive LSS, such that $ is bijective and V is 
nonperiodic. Then V does not consist of model sets, therefore V is not pure point diffractive. 

Proof. The proof is an application of Corollary 14.71 $ is bijective, thus \& admits no pairwise 
coincidence at all. All we need to show is that no ^ [ a ] ^ s a singleton set. 
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First we note that in V all types of points occur with the same frequency, if V arose from a 
bijective substitution: Since (V, <&) is bijective, every (p appears exactly once in each row of 
<E>. Since (V,<&) is admissible, every <p appears exactly once in each column of Therefore, 
the row sums of the substitution matrix S(Q) = (#&ij)i<ij< m are all equal. Let q € N be 
this row sum. Thus q is an eigenvalue of S(<&), with corresponding normalized eigenvector 
(1/m, 1/m, . . . , l/m) T . Since S(&) is primitive, it follows from the Perron-Frobenius Theorem 
that q is the unique eigenvalue of which is largest in modulus (since it is the unique one 
with a strictly positive eigenvector). It is well-known (and easy to see) that the corresponding 
(normalized) eigenvector contains the relative frequencies of the points of each type. Thus, 
all types of points have the same frequency in V . 

Assume that there is some ^q[o] such that #*$> [a] = 1. This means that all points and only 
the points in a + L' are of the same type, which occurs with relative frequency — . Therefore, 
the index [L : L'\ equals m, but there are exactly m types, and each coset L/L' contains at 
least one type. Then the only possibility is ^^> [b] = 1 for all b and V is periodic (w.r.t. the 
lattice L'). 

Under the nonperiodicity assumption of the corollary, it now follows that each ^ [a] contains 
more than one element, which completes the proof. □ 

The following theorem allows us to extend the results to nonadmissible LSS. Under a certain 
condition, we can find for a primitive LSS (V, $) an admissible one, say, (V', $'), such that V 
and V' are MLD. Here and in the sequel, V n RM d is to be read as (V\ n RB d , ...,V m n RM d ). 

Definition 4.10. Let (V , <£) be a primitive LSS and F a fundamental domain of the lattice 
supp(V). We say that (V, $) is nicely growing, if there is an R > such that, for all 
R-clusters V (1 (t + RM d ), 



In plain words: In particular, we require that the action of $ maps the set of points of 

V n RM d to the points inside a region which is considerably larger, namely, V n (QF + RM d ) 
(and maybe to some more points). And in general, the points of V inside any ball of radius 
R, centered at t, are mapped to the points in the set V n (Qt + QF + RM d ), "centered" at 
Qt (and maybe to some more points). 

Remark: All admissible substitution systems are nicely growing. If (V, $) is admissible, then 
$ maps V n RM d onto V n Q(RM d ) (and possibly more), at least for R large enough. (E.g., 
the admissible substitution of the substitution system (V, $) in Example 16 . 1 1 maps the cluster 
VHM 1 = Vn[-l,l] to the cluster Vn[-6,11], which is a superset of Vn ([0,6] + [-1, 1]) = 

V n [—1,7], thus © is fulfilled in this example for t = 0, R = 1.) Then the number of 
points in the support of the left hand side of (j2J grows approximately like R d + cR d ~ l , where 
the number of points in the support of the right hand side grows like R d \ det(Q)|. Since 
| det(Q) | > 1, Equation © is fulfilled for R large enough. This shows that all admissible 
substitutions are nicely growing. 

This is also true for the nonadmissible substitutions in Example l6.5l Essentially, all primitive 
lattice substitution systems known to the authors are nicely growing (but cf. Example I6.6|) . 
Note that the requirement of checking i?~clusters is a finite problem, since there are only 
finitely many translation classes of them. 
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Theorem 4.11. Let (V, <&) be a nicely growing primitive LSS. Then there exists (V',&), 
such that V and V' are MLD and (V',<&') is admissible. 

Note that our general assumption on V — all r-clusters in V are contained in some iteration 
$> k (v) of a single point v — ensures that V is repetitive. 

Proof. The basic idea is to construct V 1 out of V by assigning to each colored point in V a 
colored point in V 1 at the same position, where the color of the latter one depends on the 
local surrounding of the former one. This construction is well-known, see for example |KeP| 
("forcing the border") or ORS"] ("collaring"). In general, V' will use more colors than V. 
Then the admissible MFS <3?' arises in a quite natural way, using the property &(V) = V. 

Let R > such that 

supp(F) n (QF + RE d ) C supp ($(V n RS d )^ . 

Since (V, <3?) is nicely growing, we find such an R by Definition 14. 1UI Let u\, . . . , be a basis 
of the lattice L, i.e., L := (ui, . . . ,Ud)z- We introduce the fundamental domain 

d 

F = {Y^ Ajiij | < A, < 1} 
i=l 

of L. For a fixed lattice L, the lattice basis is not unique, hence there are many possibilities 
how F may look like. For simplicity, we may choose the vectors Uj such that the diameter 
diam(F) becomes minimal among all possible fundamental domains. But this is not relevant 
in the following. 

Let x S L. Then any set x + F contains exactly one element of L, namely x, and any set 
Qx + QF contains exactly | det(Q)| elements of L, which just are representatives of L/QL. 

Consider all types of i?-clusters in V, i.e., all sets 

(3) V n (x + RS d ) := (Vt n(x + RS d ), ...,V m n(x + RM d )), 

where x € L. There are finitely many translation classes of these, [C^ 1 )], [C^ 2 )], . . ., [C^ n ^], 
say. For simplicity, we represent every translation class by a cluster centered at the 
origin. Note that the are of the form (C^ , • ■ ■ , Cm ). Now we construct V' out of V in 
the following way: 

V( := {x e L | V n (x + RM d ) e [CW]}, 

and V' = (V[, . . . , V£). In plain words: If x is the center of a cluster of type in V, then x 
in V' gets color i. Since all lattice points x 6 L are covered, we get L = Ui<i<n 

V{. Clearly, 

any x £ L is the center of exactly one -R-cluster, therefore this union is disjoint. 

The MFS $' is given as follows: Let V 1 n {x} be of type j, i.e., x £ K'. Now define by 

$'(V n {x}) :=VD (Qx + QF) = ((Qx + QF) n V?, . . • , (Qx + QF) n V^). 

The sets (Qx + QF) C\V( contain finitely many points, which define the maps in <J?^- as follows. 
Let j be fixed. If (Qx + QF) n V? = {6i, 62, ... , then set := 6^ — Qx € QF n L, 
tftiv) = Qy + a%j£ (1 < < k) and $^ = {(pi, (p 2 , . . ■ , ipk}- 

We claim that (V', is an admissible primitive LSS such that V and V' are MLD. The 
maps of <£' are of the form x 1— ► Qx + a,^. By construction, the substitution given by <£' 



12 



D. FRETTLOH AND B. SING 



maps every element of V to | det(Q)| elements, which are contained in some translate of QF. 
Moreover, every is an element of QF n L, which is a fixed system of representatives of 
L/QL. Thus the definition of admissibility is fulfilled. 

Obviously, a point in V' determines an i?-cluster in V, so V is locally derivable from V' . 
By construction, V 1 is locally derivable from V, so V and V' are MLD. 

Let Vn{i} be of type i, R > 0. Since (V, <3?) is primitive, and by our general assumption (see 
the beginning of Section^, <J? fe (V D {x}) contains all types of i?-clusters for some appropriate 
k. Therefore, applying ($>') to the corresponding point V' n {x} yields all types of points in 
V. The choice of V fl {x} is arbitrary, so it can be the center of any type of -R-cluster in 
V. Altogether, this means that, for any V' n {x}, there is k € N such that (&) k (V' n {x}) 
contains all types of points that occur in V'. This is equivalent to $' being primitive. 

All we are left with to show is that <&' is well-defined, i.e., that two points of the same type are 
substituted in the same way. This follows essentially from the condition <&(V) = V: Since two 
congruent clusters in V are mapped by to congruent clusters in &(V), the corresponding 
points in V' are mapped to congruent clusters in V' . Formally, we need to show: If x, y € V- 
for some i, i.e., V' fl {x} = (V' fl {y}) + x — y, then 

(4) &(V' n {x}) = $'(V' n {y}) + Q(x - y) 

holds. So, let V' fl {x} = (V' n {y}) + x — y. By construction of V', it follows that the 
corresponding i?-clusters in V are of the same type, i.e., 

V n (x + RM d ) = {V n {y + iflB d )) + x - y. 

Since $(V") = V, we obtain 

$(yn(i| RM d )j = $ ((V n (y + iflB d )) + x - = $ n (y + iflB d )) + Q(x - y). 

Consider V D (Qx + QF + RE d ), V n (Qy + QF + i?B d ) + Q(s - y). The growing nicely 
condition implies: 

supp(Vn(Qx + QF + RM d )) c supp ($ (V n (x + iflB d ))) 

supp (V n (Qy + QF + fflB d )) + Q(x - y) C supp ($ (V n (y + + Q(z - y). 

The four sets V n (Qx + QF + RM d ), $ (V D (s + V n (Qy + QF + i?B d ) and 

<!> (V fl (y + i2B d )) occurring in the above equation are all subsets of V as multi-component 
sets, i.e., with respect to colors. We obtain 

V n (Qx + QF + c $(^n(i| iflB d )) 

Vn(Qy + QF + RM d ) + Q(x-y) C $ (V n (y + fflB d )) + Q(x - y). 

Obviously, the supports of the left hand sides are equal. Therefore, it follows from the above 
equation, that the left hand sides are equal as multi-component sets: 

V n (Qx + QF + RM d ) = Vn(Qy + QF + RM d ) + Q(x - y) 

By the construction of V' , the .R-cluster Cr = V fl (x + RM d ) in V determines the color of 
x E V . Therefore, the above equation yields 

V n (Qx + QF) = V'n (Qy + QF) + Q(x - y), 
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and by the construction of this is equivalent to @. Altogether, (V',<&') is an admissible 
primitive LSS, such that V and V' are MLD. □ 

The last result proves useful as follows. 

Theorem 4.12. Let (V,&) be a primitive LSS, and (V',&) be an admissible one, derived 
from (V, $) as in the proof of Theorem 14.111 // V' consists of model sets, then so does V . 

Proof. Let V = (Vi, . . . , V m ). By the construction of V', each set Vi is partitioned into sets 
V( , V( 2 , . . . , V- belonging to V' . Recall: An element of V' at position x got its type according 
to the i?-cluster surrounding xinV. Denote all translation classes of i?-clusters with center 
of type % by [QJ, [Q 2 ], . . . , [C ik ]. Then V = U V( 2 U • • • U V( k . 

Let V' consist of model sets. This means that each V' is a model set. To each V/ belongs a 
window set W/.. Let V( x ,V/ 2 , . . . , V( k be the sets assigned to VJ, i.e., Vi = Ui<i<S; ^7 • Let Q^. 
he the window of V' (cf-Ej). Then 

v *= U K= U ^H) = ^( U 

l<j<k l<j<k l<j<k 

Since the windows f^. have nonempty interior and compact closure, this is also true for their 
(finite) union. Thus Vi is a model set, too. □ 

5. Remarks on Complexity 

Theorem 14.51 gives an exponential upper bound for computing modular coincidences. This 
may be far away from the best upper bound. In particular, it is not known whether the 
required number of iterations is polynomial in m or not. 

If = 2 m , the A:-th power of an (m x m)-MFS $ can be computed by m — 1 "matrix 
multiplications" , to be understood as stated after Def. 12.11 If k' < k and $ fc already admits 
a modular coincidence, then also does. Therefore, it suffices to compute <E>, 3> 2 , $ 4 , $ 8 . . .. 
At first sight, this would lead to a polynomial bound for computing modular coincidences. 
However, unlike usual matrix multiplications, the costs for each multiplication is not m 3 , but 
the costs per multiplication increase with k. We do not study this question in detail here, 
but essentially, the costs grow faster than the sum of the entries of the substitution matrix 
S($ k ). 

The authors conjecture that the lowest upper bound is quadratic in m, hence not linear in 
m. Consider the following example: 

(5) 1 i— > 1 2, 2 i— > 2 3, 3^3 4, . . . , m — 1 1— > m — 1 m, mull 

If m = 2 this reads as the well-known substitution a — ► ab, b — > aa. For 2 < m < 16, the first 
modular coincidence occurs in 

q>(m-i)\ Thig 

is a hint that a general upper bound is quadratic 
in m. It will be difficult to find a worse example, since the dimension (=1) and the factor 
det(Q) (=2) are the smallest possible for admissible primitive LSS, as well as the number 
of pairwise coincidences (=1). For m £ {2,3,4,5,6}, we computed for all substitutions with 
these three properties the minimal number k, such that shows a modular coincidence (if 
there is such a k). It turned out that the substitution © yields the maximum value of k in 
these cases. 
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6. Examples of Sequences and Tilings 

6.1. Periodic substitution I. The primitive substitution given by the following MFS yields 
a periodic sequence. Here h L denotes the map hi{x) = 6x + i, i G Z. 

({ho, hi} {h ,h 3 } {h ,h 3 }\ 
<*> = {hi,h 5 } {h u h 5 } {hxM} 
\{h 2 ,h A } {h 2 ,h 4 } {h 2 ,h 4 }J 

In symbolic notation, this reads a i— > abcacb, b \— > abcacb, c i— » abcacb. This substitution on 
L = Z yields L a = 3Z, L fc = 2Z and L c = 2Z (and therefore L' = L). Note that here — as in 
some of the following examples — we use letters a, b, c instead of numbers 1, 2, 3. To obtain 
a tiling, replace the letter at position n by the interval [n, n + 1] (and label it with its letter). 
To obtain a multi-component Delone set, let V = (V a , Vf,, V c ), where V a contains all positions 
of letters a etc. 

6.2. Periodic substitution II. The trivial substitution a i— » ab, b i— > aft on L = Z yields the 
periodic solution . . . ababababababa . . . (and therefore V a = 2Z, VJ, = 2Z + 1 and V = 2Z). It 
already shows a modular coincidence relative to QL' (each < i < 3, consists of exactly 
one element). Following the construction preceding Def. 13.21 one would obtain as internal 
space H = C 2 x Z 2 (where C 2 = Z/2Z is the cyclic group of order 2, and Z 2 are the 2-adic 
integers) and the lattice L = {(t, t mod 2, t) € M x C 2 X Z 2 1 1 G Z}, but it is easy to see that 
#' = C 2 and L' = {(t, t mod 2) G M x C 2 | t G Z} suffices. 

6.3. Thue-Morse sequence. The Thue-Morse sequence is defined as the fixed point begin- 
ning by of the substitution ^ 01, 1 ^ 10 (it starts 0110100110010110100...). There 
is also an arithmetic definition of this sequence (e.g., see [Hj p. 2], |PFl 2.1.1]): For any 
integer n, denote by S 2 {n) the sum of the dyadic (binary) digits of n, i.e., S 2 (n) = ^i>o n i' 
if n = ^j>o n i^* ( n i e !})• Then the n-th symbol in the sequence is a 1 iff S 2 (n) is odd. 
Obviously, for n < 2 k , we have the following: if S 2 (n) is odd, then S 2 (n + 2 k ) is even and vice 
versa. Therefore ^fc[ n ] = I ( re + 2 fc Z) n V^} = {0, 1} for all k and n, and the Thue-Morse 
substitution admits no modular coincidence. 

6.4. Parents and children. Consider a substitution defined by a i— ► aba, b i— > 5cc and 
c i— > a6c (this substitution arises in the study of the Kolakoski sequence over the alphabet 
{2,4}, see [EH)- The corresponding MFS <3? is given by: 

/{h ,h 2 } {h }\ 
*= {^} {ho} {hi} , 
\ {h,h 2 } {h 2 }j 

with hi(x) = 3x + i. A fixed point is given by . . . ababccabcababccaba . . ., where the dot 
marks the zeroth position. It is easy to check that here we have L = Z = L'. In this case, 
as L = L' and $ is admissible, the check for modular coincidence reduces to finding a map 
g which is shows up in only one row of & k for some k G N. Therefore we get (note that 
$aa = {ho,h 2 }, $ ba = {hi}, etc.) 

m C^> aa U% b U^ ac , $[1] C$ ba U<5> cb U<S> bc , $[2] C $ aa U $ cfe U $ cc . 
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Here we use letters rather than numbers to avoid confusion, because numbers are used to 
represent cosets below. Since none of the sets <&[•] is contained in just one row of the MFS, 
itself does not show a modular coincidence. (For instance, $[0] is contained in the first and 
second row of 3>[1] in the second and last row, etc.) Nevertheless, <3? 2 shows a modular 
coincidence already, or more precisely, three of them. 

f{h%,h h 2 ,h 2 h Q ,h%} {h h 1: h h 2 } {h% h 2 h , h h 2 } \ 

$ 2 = I {h^h^h^} {h^h^h^} {h^Q, h h v h^} J 

\ {^Vi) {Mo'MoA^i^i} {huhzh^hl} J 

The map h Q o h 2 {x) = 9x + 6 is contained in just the first row, the map h 1 o h 2 (x) = 9x + 7 
just in the second row and the map h 2 o h^{x) = 9x + 5 just in the third row. 

Since L = L' , we have ^f [0] = {a,b,c}. The next step is easy: For ^[i] we just have to 
consider the maps 3>[z] (see the first indices of the matrix elements of <3? involved in to 
obtain 

* 1 [0]={o,6} ) * 1 [l] = {b,c}, * l [2) = {a,c}. 
Basically, we make use of the (unique) ternary expansion of each natural number and that 
the maps of the MFS line up with this ternary expansion. E.g., to calculate ^[7], we observe 
that 7 = 2-3 + 1, therefore ^ 2 [7] is determined by the i's of 3>y[l] with j G ^i[2], yielding 
VP 2 [7] = {b}. For completeness, we list all \E' 2 : 

V 2 [0] = {a,b}, * 2 [l]={6,c}, * 2 [2] = {a,c}, * 2 [3] = {a, b}, 

* 2 [4] = {6,c}, * 2 [5] = {c}, * 2 [6] = {a}, * 2 [7] = {6}, * 2 [8] = {a,c}. 

Here, the children of ^[i] are \P 2 [3 • i], ^ 2 [3 • i + 1] and * 2 [3 • i + 2]. Since ^ 2 [5], *2[6] and 
^2 [7] are singleton sets, we have a modular coincidence. 

Note the correspondence between & k {i] and ^jJi]- For instance, the fact that ^ 2 [5] = {c} 
corresponds to the fact that h 2 o h±(x) = 9x + 5 is contained in the last (the "c-th") row 
of $ 2 . 

6.5. Nonadmissible substitution I. Consider the one-dimensional substitution given by 
the following MFS (where hi denotes the map hi{x) = 2>x + i): 

( {h } {/ii} {h 2 } \ 
*= {hiM {h } {h u h 5 } 
\ {h 2 } {h } J 
In symbolic notation, this corresponds to the substitution 

a —> abc_Jb, b — > ba, c — ► cba- J), 
where _ stands for an empty position. A fixed point of this substitution is 
(6) ... cba bab abc bab abc bob cba bob abc . . . 

The dot marks the zeroth position. Obviously, this substitution is not admissible. Never- 
theless, a brief look at the first iterations $ fc ({0}, 0, 0) shows that this substitution is nicely 
growing. It is easy to see that this substitution is equivalent to the admissible substitution 

a — » abc, b — > bab, c — > cba. 

Both substitutions yield the same sequences, therefore the sequences are clearly MLD. Here 
we did not use the construction of Theorem 14.121 This would yield another sequence V' 
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Figure 1 . Tile substitution for the Chair Tiling (left) and for the Table Tiling 
(right). 



\ 



p r 



Figure 2. Labelling for Chair Tiling (left) and for the Table Tiling (right) 



MLD to ©, but it would use more than three letters. Namely, the construction works with 
R = 2. There are eight different types of 2-clusters (five-letter words) in ©, yielding eight 
different letters in the new sequence V . This example is easily generalized to more than one 
dimension. 



6.6. Nonadmissible substitution II. The following MFS yields another nonadmissible 
substitution: 



$ = 










{^15} \ 


{h u } 


{his} 


{hi2} 


{^14,^13} 


{his} 


{hu} 


{/ll5, /li4, /113} 





\{/ll2,M 


{hi*} 





) 



where hi(x) = 4x + i. In symbolic notation, this would look quite weird, since the action of 
$ contains large shifts to the right. E.g., a letter a at position becomes a word d- . . . -dcba 
under the action of <£, showing a large gap between the first letter d at position and the 
second letter d at position 12. This substitution is clearly not admissible. Nevertheless, there 
is a V such that $(V) = V. E.g., let C = ({-5}, 0, {-4}, 0). The iteration of $ on C yields 
a sequence with $(V) = V, namely V = UfceN & k (C). Thus (V, <&) is a primitive LSS. It is 
not easy to see whether this substitution is nicely growing. If so, it can be shown that R must 
be at least 15. (Note that the difference between the smallest and the largest translational 
part occurring in $ is 15.) 

This example can be generalized to higher dimensions, too. It can also be generalized to 
any number m > 4 of types of points such that the difference between the smallest and the 
largest translational part occurring in $ is m 2 — 1 (and we really have an LSS). In this case, 
it can be shown that R must be at least m? — 1, if the corresponding substitution system is 
nicely growing. Anyway, the substitution above on 4 letters yields the same sequences as the 
admissible substitution 

a — > dcba, b — > dbca, c — > cccb, d — > dbba, 
showing a modular coincidence in <I> 2 . 
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S p 

P Q 

s r 

r a 



P q 

s r 
p s 



s p 

q p 

r s 

r a 



q q 

p r 

p r 

s s 



Figure 3. Induced substitution for the Chair Tiling (left) and for the Table 
Tiling (right). 




{p,q,r,s} 



Figure 4. Coincidence graph for the Chair Tiling (left) and for the Table 
Tiling (right) 



6.7. The Chair Tiling and the Table Tiling. The chair tiling and the table tiling (also 
called the domino tiling) and their spectral properties are discussed in detail in |Soj . |R] 
and |BMSj . They are generated by the tile substitutions in Fig. ^ The expansion Q acts as 
scaling by a factor 2, and the rotated versions of the shown tiles are substituted analogously. 
Both tile substitutions can be described by substitutions over the alphabet {p, q, r, s} on 
L = Z 2 , where one uses the labelling of Fig. |^1 

These labellings yield the (induced) substitutions in Fig. 03 Therefore, we obtain the coin- 
cidence graphs shown in Fig. 0] (observing that L = V = Z 2 ). Note that the substitution 
for the table tiling admits no pairwise coincidence, therefore Corollary 14.71 yields that the 
induced LSS does not consist of model sets. The substitution is also a bijective substitution, 
so Corollary 14.91 applies, too. 

These graphs show that the chair tiling is pure point diffractive (compare to O Theorem 
8.1], also see |So| Example 7.1]), while the table tiling is not (compare to _R, Theorem 8.2], 
also see |So| Example 7.3]). We note that the internal space for the chair tiling is Z 2 (which 
lines up with the inflation by a factor of 2, i.e., Z| = hm^ I? /2 k ■ Z 2 ). 
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Figure 5. The substitution for the Semi-Detached House Tiling (left) and a 
part of the tiling, which (hopefully) illustrates its name. 

6.8. Rl R2 Tilings. In |Fra1| . a family of substitutions is considered which fulfil two 
conditions, namely, (Rl) and (R2). In our notation, these read as follows. Let (V,<&) be an 
admissible primitive LSS supported on L = V = (V\, . . . , Vim), [L : L'\ = m. 

(Rl) ^ol *] = {*) * + m } (1 — * — m )> wriere a i runs through L/L' . 
(R2) n ^.i+m = for all 

The first condition tells us that no ^q[o\ is a singleton set: each one contains two elements 
i,i + m. Therefore, the second condition tells us that there is no pairwise coincidence in any 
^ofa]. Thus, it follows from Corollary 14.71 that none of the multi-component Delone sets (or 
tilings) induced by these substitutions consist of model sets. 

6.9. The Semi-Detached House Tiling. The tile substitution in Fig. |S] (with expansion 
factor 2) appears as unsolved problem in [Fl Abbildung 17]. As above, rotated and reflected 
versions of the tiles are substituted analogously to the shown tiles. 

Each of the two semi- houses (i.e., the big and the small one, see tile substitution in Fig. EJ) 
arises in 4 different orientations and its reflected versions. To derive a substitution, we use 
the labelling in Fig. H3 

We then obtain the substitution for "squares" in Fig. [71 p. HOI which is again a substitution 
on L = 1? by a factor of 2. 

Instead of using all 56 squares which occur in this substitution, we form 20 classes which are 
already present in Fig. [7J (i.e., squares which yield the same substitute are in the same class). 
We now look at the substitution induced on these classes, e.g., 

i e _ e i 

a — > _ a — > _ . 

a a a a 

Obviously, if we do not have a modular coincidence for this induced substitution, then there 
is no coincidence for the original substitution either. The coincidence graph for the induced 
substitution is shown in Fig. El p. As suggested by Fig. UJ L' forms a checkerboard with 
respect to L, i.e., [L : L'] = 2. 
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Figure 6. Labelling of the Semi-Detached House Tiling 



As we do not have a modular coincidence for this substitution (no vertex in the coincidence 
graph represents a singleton), the semi-detached house tiling does not arise from a model set 
and is not pure point diffractive. This non-trivial example shows that analyzing LSS and 
tilings can often be quite elaborate. 

Appendix A. Modular Coincidence vs. Dekking Coincidence 

An admissible LSS in one dimension (i.e., on L = Z) is known as a primitive substitution 
of constant length q, see [DJ [^J IPFj . Here the mappings of the MFS <E> are of the form 
x i— »• q ■ x + a- k where a- k G {0, . . . , q — 1}, i.e., each letter is substituted by exactly q letters. 

The first step in analyzing such a substitution is to determine L'. Since L = Z, we have 
L' = r ■ Z and r E N is to be determined. To this end, we define, for a sequence u = 
. . . u_ 2 u_ 1 u u 1 u 2 ■ ■ ., the sets S k = {j € Z | u- +k = u k } and numbers g k = gcd S k . It is easy 
to see that g k = g k , if u k = u k , , so that we only have to determine the g k once for each letter 
of the alphabet. Furthermore, we have L- = g k ■ Z if u k = i (in other words: if is of type i) 
and, by L' = L 1 + . . . + L m (m is the cardinality of the alphabet), also r = gcd{g k \ k G Z}. 

Let us introduce the height h of a substitution of constant length q (see [DJ Definition II. 8] 
and Q, Definition VI. 1.]): h = max{n > 1 | gcd{n,q} = 1 and n divides g }. We have the 
following property for every k € Z ([Dl Remark II.9(i)] and [^J VI. 2. 1.1.]): 

{n > 1 | gcd{n, q} = 1 and n divides g } = {n > 1 | gcd{n, q} = 1 and n divides g k } 

From this, it follows directly that 1 < h < r < m and h divides r, but is relatively prime to 
q. If r = m, then the sequence under consideration is periodic (compare to the argument in 
the proof of Corollary I4.9j) . It can happen that h < r. 

Example: The paperfolding sequence can be generated by the substitution a i — > a6, 6 i — ^ cb, 
c i— > ad, d i — ► cd and a subsequent projection {a, b} *— > 0, {c, d} i— ► 1 on a two letter alphabet 
{0, 1}. We consider the original substitution on {a, b, c, d} for which we have: q = 2, g k = 2 
for all /c G Z, r = 2 and h = 1. 
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Figure 7. Substitution on squares, derived from the substitution in Fig. [21 



For a substitution a of constant length with h > 1, it is always possible to introduce a 
substitution r\ on a (usually bigger) alphabet of the same constant length with height 1 such 
that the sequence generated by a can also be obtained from the sequence generated by r/, 
see [01 Lemma 11.17 &: Definition 11.20]. We illustrate the general method by an example, 
also compare this construction to Theorem 14. Ill 

Example: Consider the substitution (see [DJ Example following Definition 11.20]) 
h- > 010, 1 h- > 102 and 2 i— ► 201 of height 2. This substitution generates a sequence u = 
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Figure 8. The coincidence graph of the derived tiling of the Semi-Detached 
House Tiling. 

. . . 010102010102 .... Now we take as the new alphabet all words in {u n , h . . . u n . h+h _ 1 \n € Z}, 
in this case just a = 01 and b = 02. Then the substitution n of height 1 is given by a 1— > aab 
and b 1— > aba (since a = 01 1— > 010102 = aab). 

Let n be a substitution of constant length q and height h = 1. We say ?? admits a coincidence 
(or Dekking coincidence) if there exist /c and < j < q k such that the j-th symbol of r/ fc (i) 
is the same for all letters i ([QJ Definition VI. 3], see also [0 Defnition III. 1] ) . In the above 
notation, 77 admits a coincidence if n k (i)j is the same for all i. A substitution a of constant 
length q and height h > 1 is said to admit a coincidence if its associated substitution of height 
h = 1 admits one. 
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Reformulating [HJ Theorem 11.13 & Theorem III. 7] (compare to [QJ Theorem VI. 13 & The- 
orem VI.24] and [PFj Theorem 7.3.1 & Theorem 7.3.6]) we obtain: 

Lemma A.l. Let a be a substitution of constant length q and height h. Let u be the sequence 
generated by a. Then u is a regular model set iff a admits a Dekking coincidence. The internal 
space H is given by H = r L q x C h , where TL q denotes the product over the distinct primes p 
dividing q of the p-adic integers 7L, p and C h is the cyclic group of order h. □ 

We compare this result with the ones involving the modular coincidence (Definition 13.21 and 
Theorem 13. 3[) : The substitution has a modular coincidence iff it is a regular model set with in- 
ternal space C r x lim^fc Z/g fc Z. Obviously, lim^fc Z/g fc Z ~ Z (e.g., observe that lim<_fc Z/4 fc Z 
is topologically isomorphic to Z 2 , wherefore it is enough to consider just the distinct primes 
p). We can write r = h ■ q, where gcd{r, q} = q, and we observe that C r x Z is topologically 
isomorphic to C h x 7L q (note that C~ x 7L q and "L q are topologically isomorphic). So both 
methods actually give - up to isomorphism - the same internal space. 

In a sense, the height h tells us whether L' = r-Z lines up with the length q of the substitution. 
If h = 1 but r > 1, then we do not have to know the sets 3> [-] explicitly for the Dekking 
coincidence, the substitution rule itself already gives the right partition of the alphabet. 

Example continued: The underlying substitution of the paperfolding sequence (a i— > ab, 
b i— > cb, c i— ► ad, d i— ► cd) has height h = 1 but r = 2. We observe that there is either an 
a or a c on the zeroth place of the substitute, while on the first we have either b or d. But 
we also have ^q[0] = {a,c} and *q[1] = {b, d}, therefore the substitution lines up with this 
partition. We note that this substitution has a modular coincidence relative to 2 2 Z (after 
two substitutions the first letter is always an a, a i— > ab i— > abcb, b i— > c6 i— > adc& etc.), and 
therefore is pure point diffractive. 



Appendix B. Subdivision Graph for Overlaps 

(2) 

Instead of the coincidence graph G^ v ^, we can also define a pair coincidence graph G) v ^y. 

(2) 

(V$) ^ S a directed (multi-)graph, where each set of pairs of elements of ^ [ a ] f° r some a 
(i.e., both elements are from w ith the same a) is represented by a vertex. We call pairs 

of the same element {i,i} = {i} a coincidence in the pair coincidence graph. We draw an 
edge from v 1 to v 2 iff v 2 is a "child" of v 1 (Here we mean by an "child": If v 1 = {s 1 , . . . , s^\ 
then v 2 = {^(s 1 ) z , . . . , <3?(s^) 2 } for some z.) The pair coincidence graph of the table tiling is 
shown in Fig. 

Obviously, if the coincidence graph G^ v ^ has a singleton set (i.e., we have a modular coinci- 

(2) 

dence), then from each vertex of G^y ^ there is a path leading to a coincidence. Conversely, 

(2) 

since the alphabet is finite, if from each vertex of G,y ^ there is a path leading to a coin- 
cidence, then Gfy ^ also has a singleton set (compare this to Corollary 14 .7j) . Therefore, we 
can also use the pair coincidence graph to decide the question about modular coincidences 
(the pair coincidence graph of the table tiling is not connected, the coincidences forming one 
connected component; the table tiling has no modular coincidences). 

On the tiling level, the pair coincidence graph is called subdivision graph for overlaps, see |So| 
p. 721] and |LMS2| Appendix 2]. The construction of this subdivision graph follows closely, 
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o 

Figure 9. Pair coincidence graph for the Table Tiling 




Figure 10. Subdivision graph for overlaps for the Table Tiling 

with the necessary changes, the construction of the pair coincidence graph. For the pair 
coincidence graph we can think of placing two multi-component Delone sets on top of each 
other (or, placing a translated (by a vector of L') copy of the multi-component set on top) 
and observing what happens under substitution. Likewise, we can think of placing two tilings 
on top of each other and than we observe, what happens to the overlaps of the tiles under 
substitution. As an illustration we give the subdivision graph for overlaps (overlaps are 
denoted by filled areas) for the table tiling in Fig. ^] (where we identify rotated versions of 
the configurations), see also |Sol Fig. 7.4 (d)]. This graph should be compared with the pair 
coincidence graph in Fig. (basically, rotation identifies all the singletons, {p, r} and {q,s}, 
and the four sets {r, s}, {q, r}, {p, s} and {p, q}). 

Appendix C. Connections to Automatic Sequences 



In [H], von Haeseler is interested in the question, under which conditions a sequence is auto- 
matic and/or generated by a (not necessarily primitive) substitution of constant length (and 
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in the higher dimensional case even whether there exists an LSS). To this end, he introduces 
the kernel of a sequence u = uqu±U2 ■ ■ ■ jE] Definition 2.2.18.], which (in one dimension) one 
can think of as the set {u a u q k +a u 2q k +a . . . \ k G No, a G {0, . . . , q k — 1}}, i.e., the set of subse- 
quence on q kr L + a (Q k L + a) for all k and a. The finiteness (and cardinality) of this set is in 
question, and can be answered (see [Hj Lemma 2.2.22.]) by the substitution graph [HJ Defini- 
tion 2.2.20.]: The substitution graph is the directed, labelled graph with base point (1, . . . , m) 
(the whole alphabet) and the vertices are given by all ordered m-tuples (l) z /, . . . , <E> fc (m) 2 ,) 
(where z ! G L/Q k L), which can be obtained by some power k of the substitution. We draw 
an edge with label z from v l = (s 1 , . . . , s m ) to v 2 = (^(s 1 ) z , . . . , <fr(s m ) z ) (z G L/QL). 

Obviously, in the case L = L', the coincidence graph can be obtained from the substitution 
graph: m-tuples which are equal as sets are identified. And Lemma 14.31 tells us that with 
this identification there is no problem with the edges. Furthermore, an m — tuple (i, . . . , i) 
indicates a coincidence. As an example, we give the substitution (sec H, Example 1. following 
Definition 2.2.20.]) and coincidence graphs of the Thue-Morse sequence (o i— > ab, b i— > ba). 




Figure 1 1 . Coincidence graph (left) and Substitution graph (right) of the 
Thue-Morse substitution 
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